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Abstract 

In this paper, we answer a question raised by Len Pitaevskii and prove that the 
ground state of the Gross Pitaevskii energy describing a Bose Einstein condensate 
at low rotation does not have vortices in the low density region. Therefore, the first 
ground state with vortices has its vortices in the bulk. This is obtained by proving that 
for small rotational velocities, the ground state is multiple of the ground state with 
zero rotation. We rely on sharp bounds of the decay of the wave function combined 
with weighted jacobian estimates. 

1 Introduction 

Among the many experiments on Bose Einstein condensates, one consists in rotating the trap 
holding the atoms in order to observe a superfmid behaviour: the appearance of quantized 
vortices [1, 23, 19, 20, 2]. This takes place for sufficiently large rotational velocities. On the 
contrary, at low rotation, no vortex is detected in the bulk of the condensate. The system can 
be described by a complex valued wave function minimizing a Gross Pitaevskii type energy. 
A vortex corresponds to zeroes of the wave function with phase around it. The density 
of the condensate is significant in a region which is either a disk or an annulus, and gets 
exponentially small outside this domain. Vortices are experimentally visible in the bulk of 
the condensate. A question raised by Len Pitaevskii is whether for small rotational velocity, 
when there are no vortices in the bulk, vortices could exist in the low density region. For 
very large rotational velocities, when bulk vortices are arranged on a triangular lattice, it 
has been shown [4] that the vortex distribution extends to infinity. Therefore, in this case, 
there are many vortices in the low density region. It is then very natural to wonder whether 
vortices first appear in the bulk or at infinity. It is experimentally and numerically difficult 
to observe a vortex, which is a zero, in a low density region. Mathematically this could not 
be achieved through energy estimates or expansion since the contribution of a vortex in a 
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low density region is very small. In this paper, we introduce new ideas to answer Pitaevskii's 
question and prove that at low velocity, there are indeed no vortices in the condensate, even 
in the low density region. Therefore, the first ground state with vortices has its vortices in 
the bulk. 

Since a condensate is a trapped object, the geometry of the trap plays a role. An 
important special case is a radial harmonic trapping potential V(r) = r 2 . The space can then 
be split into two regions, a region of the form V = {A > V(r)} (for a suitable constant A ), 
where the wave function is significant and the condensate is mainly located, and a region 
R 2 \ V where the modulus of the wave function is exponentially small [2]. In this latter 
region, it is very difficult to determine mathematically the contribution of a vortex to the 
energy. Ignat and Millot [11, 12] have determined the critical rotational velocity fl c for the 
nucleation of the first vortex inside T>. This theorem does not describe the behaviour in 
R 2 \ T>. A natural question is whether for Q < Q c , the minimizer of the energy has zeroes in 
this region, whether there is a smaller critical velocity than Q c where the minimizer is unique 
and vortex free. At very high velocity, it has been proved in [4] that vortices exist up to 
infinity so it seems reasonable that at smaller velocity, vortices may exist in the exponentially 
small region, far away from the bulk and could arrange themselves on disks or arrays close 
to infinity. In fact, we prove that this is not the case before f2 c , namely that the minimizer 
is unique and does not vanish. It means that for a large range of rotational velocities Q, the 
minimizer is given by the same function. 

We consider here a two-dimensional setting and define the energy for the complex- valued 
wave function u, such that / R2 \u\ 2 = 1, as 

where Q is the angular velocity, x = (xi,Xi), x L = (—X2,Xi), e > is a small parameter, 
V(x) is the trapping potential and (iu, Vu) = iuVu* — iu*Vu. The class of potentials 
includes the model case V = x\ + x\. Then, the critical angular velocity for nucleation of 
vortices is of order | loge| (see [11]). An upper bound on the rotational velocity is given by 
Q < 1/e when the confinement breaks down. The condensate is mostly concentrated in the 
region 

V := {x E R 2 : V < A } (1.2) 

where Ao is chosen so that 

(A - V(x)) + dx = 1. (1.3) 

We refer to [2] for more details on how this is derived from the physical experiments. 

In recent experiments in which a laser beam is superimposed upon the magnetic trap 
holding the atoms, the trapping potential V(x) is of a different type [21, 23, 24, 5]: 

V(r) =r 2 + V e~ r2/wo . (1.4) 

When the gaussian is expanded around the origin, this leads to a harmonic plus quartic 
potential [23] 

V(r) = (1 - b)r 2 + k -r\ (1.5) 
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If b is small (b < 1 + (3/c 2 /4) 1//3 ), the domain V given by (1.2) is a disc, while if b > 
1 + (3/c 2 /4) 1//3 , it is an annulus. According to the values of Vq and wo in the case of (1.4), 
the domain V can also be a disk or an annulus. 

In this paper, we consider potentials V including r 2 and of the type (1.4) or (1.5) when 
the bulk V is a disk. In the case where V is a disk, the potential V is not necessarily required 
to be increasing. 

1.1 Assumptions 

Throughout this paper, we make the following assumptions about the potential V . First, 

V is nonnegative and radial, V G C 1 , (1.6) 

and 

there exists c > 0,p > 2 such that — r p < V(r) < c r p if r > c . (1.7) 

c 

This assumption is easily seen to imply that E £ is bounded below for \Q\ <C 1 and that the 
angular momentum term x 1 - ■ (iu, Vw) is integrable as long as u has finite energy. We will 
also use (1.7) to obtain decay estimates that justify for example the integration by parts 
leading to a decoupling of the energy. We fix Ao G R such that (1.2)-(1.3) hold. Such a 
A exists due to the growth of V. We further assume that the bulk V is a disk and not an 
annulus, that is V is such that 

V = B R {0) for some R > (1.8) 

and that there exist 5 > and a C 1 function R : (— 25 , 25 ) — > R also denoted R$ = R(S), 
such t licit 

Ro = R, {x:V(x) < \ + 5} = B Rg (0) 

(1.9) 

and < ^ < dR/d\ < C for some constant C 

where B r (y) denotes the open ball of radius r about y. This implies that Ao — V is bounded 
away from in the interior of V; in physical terms, this assumption rules out the case of 
annular bulks and "giant vortices" at low angular velocities. We remark that the assumption 
above implies that if \x\ G (R-s, Rs) and < 5 < 5q then dist(x, &D) = 0(5). 
We point out that assumptions (1.7) and (1.9) imply that 

there exists c x > such that V(r) - A > cx(r 2 - R 2 ) for all r > R. (1.10) 

Our assumptions include indeed potentials like r 2 or (1.5) for a disk case, and do not 
require V to be increasing. 

1.2 Main result 

Our main result is 
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Theorem 1.1. Assume that u £ minimizes E £ (-) with rotation Q, and let r\ £ denote the 
minimizer of E £ {-) for Q = 0. There exists eo,oJo,u\ > such that if < e < Sq and 
Q < u \ loge\ — u> 1 log | log e| then u £ = e %a r\ £ in R 2 for some constant a. 

In the pure quadratic case V = r 2 , Ignat and Millot [11, 12] have shown that the bulk 
of the condensate (that is any domain contained in V) is vortex- free for \Vt\ < u \ log e| — 
u>i log | logs], for some oj\ > and the same constant u that we find in Theorem 1.1 . They 
have no information on what happens in R 2 \ T>. Our theorem proves that vortices do not 
lie in R 2 \ T>. They have also shown that there exists 5 > such that the ground state has 
at least one vortex in the bulk if Q > uj \ loge\ + 5 log | log e | . In this sense, our estimate 
\fl\ < wo I logs | — ^1 log I loge| captures the sharp leading-order term, and the correct scaling 
of the next-order term, of the critical velocity for vortex formation. We point out that our 
arguments also deal with more general potentials. The arguments used in [11] to prove 
the existence of interior vortices for rotations greater than co>o|loge| + <5 log | log £ | should 
extend with few changes to the more general potentials considered here, using results about 
auxiliary functions that we establish in Section 3 in place of parallel results from [11]. Thus 
the constant u should also be sharp for these more general potentials. 

We split the proof into two independent results. The first main result of this paper 
asserts roughly speaking that symmetry breaking occurs first in the interior of T>: if Vt is 
small enough that there are no vortices in V, then there are no vortices anywhere. 

Theorem 1.2. Assume that u £ minimizes E £ (-) with rotation Q, and let i] £ denote the 
minimizer of E £ {-) for Q = 0. Assume also that Q < C\ log e| for some C . 

There exists e > such that if < e < e and Q is subcritical in the sense that 

K\ > \r} £ m £>i := {x e V : dist(x, dV) > | loge|" 3/2 } (1.11) 

then u £ = e ta r] £ in R 2 for some constant a. 

Our second main theorem gives an estimate for the critical value of Q. The statement of 
the theorem refers to an auxiliary function f : let 

a(x) = A - V(x), 

roc r q if r > i? 

Vo:=Va+, Ur) = J r f (r) = | &(r)/|J g (r) if r < R . (1-12) 

Theorem 1.3. Let ujq = . There exists ui 1 > and S\ > such that if \Q\ < 

1 1 JO 1 1 oo 

Wo I log ^| — oo\ log I log ^| and < e < E\, then Q is subcritical in the sense of (1.11), and the 
conclusion of Theorem 1.2 thus holds. 

In our proof of Theorem 1.3, as in estimates of the critical rotation in works such as 
[11] and [3], a main point is to obtain sharp energy lower bounds. In all earlier works that 
we know of, this is done using the vortex ball construction originally introduced by [13] 
and [22]. In our proof of Theorem 1.3, we avoid any explicit 1 mention of vortex balls by 

1 However, the proof of Lemma 4.1, see Lemma 8 in [14], ultimately relies on a vortex ball construction 
appearing in [15]. 
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instead appealing to a result from [14], stated here as Lemma 4.1. This makes our argument 
considerably shorter than those in [3, 11] and other references. 

We point out that the results of [11, 12] do not directly imply that Theorem 1.3 holds in 
the case V = r 2 , although it is possible that this conclusion can be extracted with relatively 
little effort from arguments in these references. 



1.3 Main ideas of the proof 

The energy minimizers with Q = provide real solutions to the Euler-Lagrange equations: 
when f2 = 0, E £ (rj) = G £ (rj), where 

Gs(v) = / R2 {^| 2 + ^M 4 + ^V(x)\ V \ 2 ] dx. (1.13) 

Our main goal consists in proving that up to the critical velocity of nucleation of bulk 
vortices, the minimizer of E £ with velocity Q is in fact equal to i] £ . 

The minimizer r] £ of G £ under the L 2 constraint of norm 1, is (up to a complex multiplier 
of modulus one) the unique positive solution of 



-A V£ + ^Ve(V(x) + rj 2 ) = ^X £ r) £ (1.14) 



where \\ £ is the Lagrange multiplier, which is also necessarily unique. Moreover, A £ — > Ao, 
and rj 2 converges to a + in L 2 (T>) and uniformly on any compact set of T>. We will need some 
estimates on the decay of i] £ at infinity that we prove in section 2. 

By a remarkable identity (see Lassoued & Mironescu [17]), for any u, the energy E £ for 
any Q splits into two parts, the energy G £ (r] £ ) of the density profile and a reduced energy of 
the complex phase v = u/i] £ : 

E £ (u) = G £ (rj £ ) + F e (v) (1.15) 

where F £ (v) = j^^\Vv\ 2 + ^(\v\ 2 - l) 2 - r] 2 ^ ■ (iv,Vv)^j dx. (1.16) 

In particular the potential V(x) only appears in G £ . We will recall the proof of (1.15), as 
well as that of (1.18) below, in Section 3. This kind of splitting of the energy is by now 
standard in the rigorous analyses of functionals such as E £ . 
Next, define 

s V 2 (s) ds, (1.17) 
so that V" 1 ^ = x ± 7] 2 . An integration by parts yields 

Fs(v) = I'M flV^I 2 - ^Jv) + ^(H 2 - I) 2 ) dx (1.18) 



/R2 [ 2 \ T] 2 J As 2 

where Jv = |V x (iv , Vi> ) = (iv Xl ,v X2 ) is the Jacobian. 

We recall that the function f £ := C, £ /r] 2 appearing in F £ is important since it is well 
known that vortices in the interior of V first appear near where this function attains a local 
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maximum [2, 3, 11, 12]; its importance is also clear from (1.18), since it controls the relative 
strength of the positive and negative contributions to F £ . The proofs of Theorems 1.2 and 
1.3 rest on new bounds for f £ in R 2 \ V and near dV, which in turn rely on decay estimates 
for r] £ . In particular, we show in Lemma 2.4 that f £ < Ce 2 ^ in R 2 \ T>. 

The other part of the proof consists essentially of bounds of 2Q J i] 2 f £ Jv by the positive 
terms in F £ . Away from the bulk, we use our estimates of f £ to find that 2Qf £ Jv is bounded 
pointwise by ||Vt> | 2 . In the bulk, where r/ 2 is not too small, we have 

\v 2 £ Wv\ 2 + ^- 2 (\v\ 2 - I) 2 > , 2 [i|V,| 2 + ±(\v\ 2 _ lf] 

for some e such that | \oge\ = | loge|(l + o(l)). We obtain the desired bounds by combining 
this with a weighted Jacobian estimate mentioned above, Lemma 4.1, which directly implies 
that 

2Vt J xVsfeJv < n ^M^^j J xr f £ [^\W\ 2 + ^(M 2 - I) 2 ] + small error terms 

where x is a cutoff function supported in the bulk. Note that the leading-order critical 
rotation oo is such that f2 ( 2 ^Hp ?s Q/ujq\ log £ | . The proof of Theorem 1.3 is completed 
by assembling these ingredients and controlling error terms. The proof of Theorem 1.2 relies 
on an additional ingredient, which is that if |i>| > \ in an open set U, then Jv is extremely 
close in U to ^(j^) = 0. Theorem 1.1 follows immediately from combining Theorems 1.2 
and 1.3. 

An interesting open problem is to see to what extent this analysis continues to hold if the 
assumption of radial symmetry is dropped. In our arguments, this symmetry is used heavily 
in our analysis of the behavior of f £ away from the bulk, and near the boundary of the bulk. 

We briefly remark on the assumption (1.7) of quadratic growth. Our proofs show that the 
absence of vortices in the low density region is a consequence of the fact that the auxiliary 
function f £ = C, £ /i] 2 is very small in R 2 \ T>. The proof of this fact (see see Lemma 2.4) can 
be modified to show that if for example (1.7) holds with p < 2, then f £ (r) > Cer l ~' p l 2 — > oo 
as r — > oo. However, in this situation E £ is unbounded below for any Q ^ 0. This reflects 
the fact that a subquadratic trapping potential is not strong enough to contain a rotated 
condensate. 

2 Properties of auxiliary functions 

In this section we study the real-valued minimizer i] £ and the auxiliary functions £ e and 
fe = is hi defined as 

Ur) = / srf £ {s) ds, f e (r) = £ £ (r)A? 2 (r). (2.19) 

J r 
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Theorem 2.1. Assume that V satisfies (1.6), (1-9). Then for every e > 0, there exists a 
unique positive minimizer r\ £ of G £ in 

U := {u e H\R 2 ) : / \u\ 2 V(x) < oo, / \u\ 2 = 1}. 

JR. 2 JR 2 

Every minimizer of G £ in H has the form e ta i] £ , for a constant. Moreover, i] £ is a radial 
smooth positive function and satisfies (1.14) 

\X £ -X \<Ce\loge\ 1/2 (2.20) 

where Ao is defined by (1.3). Finally, recall the notations Rs from (1.9) and a = Ao — V, the 
following estimates are satisfied: 

Ve(r) < Ce^e"" 1 "^-^ in~R 2 \V (2.21) 

\r] £ -VaJ\ < C£ 1/3 v^+ inB R (2.22) 

||V77 £ ||l-(r 2) < Ce-\ (2.23) 

rf e (r) < for all r E (R_s , R$ ) (2.24) 



|^s( r )| < ~ r leij r )\Jv( r ) for all sufficiently large r (2.25) 



Certain parts of the proof follow quite closely arguments given in [3] and in the pure 
quadratic case in [11]. Note that some arguments in [11] rely strongly on the special shape 
of the potential and cannot be generalized to other functions. Since V is not necessarily 
increasing, we have property (2.24) only in the neighborhood of &D. 

Proof. Step 1. Existence of minimizer s: This follows from standard arguments once we no- 
tice that J R 2 \u n \ 2 Vdx < C is uniformly bounded for any sequence (u n ) minimizing G £ , and 
the set of functions in T-L satisfying such a uniform bound is precompact with respect to 
weak convergence in iJ 1 (R 2 ). This last fact is proved by straightforward and well-known 
arguments, such as are explained in the proof in [11], Lemma 2.1, for V quadratic, the point 
being that the bound on / \u\ 2 V prevents mass escaping to oo. Standard theory then implies 
that any minimizer is smooth. If 7] is any minimizer, then \q\ is as well, since (j(|C|) < ^(C) 
for all (. The strong maximum principle then implies that \r]\ (and hence 77) never vanishes, 
and since G{rj) < G{\rj\), it is easy to see that r)/\r)\ = e ia for some constant a. We henceforth 
let r) £ denote a fixed positive minimizer. 

Step 2: uniqueness of r\ e . Multiplying (1.14) by i] £ and integrating by parts we find that 
ji £ is positive. Suppose that there are two couples (770,^0) and (771, /ii) satisfying (1.14) such 
that 1 1 770 1 1 l 2 = 1 = II^iIIl 2 and /i > /ii, and define w — ^. This function verify 

/ r] 2 (w — l) 2 dx — 2 (qf — r/orji) dx = 2 r] 2 w(w — 1) dx 
>/Pt^ i/R,^ i/rt^ 



A. Aftalion, R.L.Jerrard, J. Royo-Letelier 



8 



-V • (f| 2 Vw) + -^r]lw(w 2 - 1) = (/I! - fi )Vo w - 



and 

1 

7 2 

Multiplying the second equality by (w — 1), integrating by parts and then using the first 
equality we find 



/ R2 {vl\V(w - 1)| 2 + ^vfai™ - 1)> + 1) + - l*M(w - l) 2 } rfx = 0. 

The integration by parts is justified in view of (2.21), (2.25), which apply to both 770 and rji, 
and the proofs of which do not rely on the uniqueness of the minimizer. Hence w = 1 and 
Ho = Hi. 

Step 3: estimate of \ e — Ao- We next note, following standard arguments, that G £ can 
be rewritten 

GM = I {i|V,|' + ^(,* - + ^ V} «fc + ± (A - i /<„♦)>) 

if H77H2 = 1. Let G\(rj) denote the first integral above. We claim that 

Gj(ife)<C|loge|. 

Since rj e is a minimizer, to prove this it suffices to construct a competitor for which G\ is 
suitably small. To do this, define 

Note that 

Using this and explicit calculations such as those in [14], Lemma 12, the claim is easily 
verified. We now multiply (1.14) by r] £ , integrate by parts and rewrite, recalling the L 2 
constraint, to find that 

j 2 {K - A ) = / |V?? £ | 2 + ^(rj 2 £ - (V - X ))rj 2 dx (2.26) 

= J \^Vs\ 2 + -^(Ve ~ a+ + a ~)Ve dx 

= J |V»fe| 2 + J 2 [arrfc + (Ve ~ a + ) 2 + (i£ - « + )« + ] (2.27) 
< 4Gl(rk) + ^\\Ve ~ a + \\ L 4a + \\L* < C[G\{ru) + -^G]^)]. 
Thus we have proved (2.20). 
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Step 4 ■ estimates ofrj £ . 
We claim that 

7] 2 £ < max(A £ - V) =: A (2.28) 

To see this define w — ^(rj £ — y/A). We have that rj £ e Lf oc , so after (1.14) w, Aw £ Lj oc . 
Kato's inequality gives Aw + > sgn + {w)Aw. Using (1.14) again we find 

£ £ 

Hence we have -Au> + + (u> + ) 3 < in £>'(R 2 ) and w G L 3 oc , so using Lemma 2 in [9], 
w + = 0. 

We remark that the properties of the potential V at the boundary (1.9) implies that the 
maximum of A e — V is attained at an interior point x of V such that dist(x , &D) > c5 . 

The minimizer being a solution of (1.14) in L°°, by elliptic regularity we derive that it is 
a smooth function. 

Proof of (2.21). We construct a supersolution of (1.14) of the form 

' ^/Ao - V(x) + 85 if \x\ < R_ 5 
fj( x ) : = J X °- 5 -J s {x) + 3V6 if < \x\ < R s 

7 e if -R<$ < |x| 

where < 5 < 5 is small parameter that will be determined later and 7, a are chosen such 
that fj E C^R 2 ), i.e., 

7 = exp ( a R s ' J and a = 16 



|W(f?a 



A straightforward computation shows that for 5 = Ce 1 ^ 3 , r\ is a supersolution of (1.14) and 
we also have 

a = G(£ 1 / 3 ) and 7 = (e 1/6 e^ 1/3 ^) . 
Moreover, with this choice of 5, ff > \ £ — V for every \x\ < Rs, so using (2.28) 

r) 2 e {x Q ) <A = \ £ - V(x ) < fj(x ). 
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Because rj £ and 77 are going to zero at infinity, if the function r\ £ — 77 is positive somewhere 
in (r , oo), for r := |xo|, then it attains a positive maximum at f G (r , oo), i.e. rf e {f) = fj'(f) 
and r)"(f) < fj"(f). Given the structure of (1.14) and because fj is a supersolution and r\ £ a 
solution, if V(f) — \ £ > we would have that rj £ (f) < 77(f). In another hand, if V(f) — \ £ < 
then we would have 77(f) < \J\ — V(f), which for e small enough, contradicts the definition 
of 77. Hence 

Ve(r) < fj(r) in (r , 00) . 

Proof of (2.22). Using assumption (1.9), by exactly following [3], one finds that |?7 e — 
< Ce 1 / 3 ^, for a £ := X £ - V = a + A e - A . In view of (2.20), this implies (2.22). 

Proof of (2.23). For x G R 2 define r/(y) = r] £ (e(y — x)) in B 2 l{x). This function satisfies 

Af) = 77 (V(e(y - x)) + fj 2 - \ £ ) =: h £ 

After estimates (2.21) and (2.22) \h e \ < C, so using a Holder estimate for the first derivative 
of fj (see Theorem 8.32 in [10]) we have that || V l fj ] \i l o °{B L (x)) < C for a constant C indepen- 
dent of x and hence the result. 

Step 4 : Proof of (2.24). 

We denote L the elliptic operator obtained by linearizing equation (1.14) 

L:=-A + \(V(x)+3r ] 2 £ -X £ ), 



e 

and Xj, j — 1, 2, ... , its eigenvalues in R 2 . 

Let fi be the first Dirichlet eigenvalue of L in the half space Q = {xi > 0} and ip the 
corresponding eigenfunction (which exists because of the compact embedding of % in L?). 
Since V and % are radial, is clear that the odd extension of ip to R 2 is a eigenfunction for 
L in R 2 with corresponding eigenvalue jj, = \j. Note that j > 2 because the odd extension 
change sign in R 2 . 

We have that Li] £ = 2i] £ > and 7] £ > 0. Using the maximun principle due to Berestycki, 
Nirenberg and Varadhan [6], this implies that the first eigenvalue of L is positive. We will 
prove that if (2.24) does not hold, then /i < 0, which contradicts the fact that Ai > 0. 

Assume that rj' £ (r) > at some r G (R-s , Rs )- Then there exists a < r < (3 such that 
rj' £ (a) = rj' e (/3) = and i]' e > in (a, (3). If a < R-25 , then r\ £ is increasing on (R_ 2 s , R-s ), 
so that rj £ (R_28 {) ) < Ve(R-8 )- This is impossible for all sufficiently small £, since r\ £ — > y/a^ 
uniformly for r < i?_ e i/ 3 , by (2.22), and a + (i?_ 2< 5 ) > a + (i?_5 () ). Thus a > R-25 - The same 
argument, but using (2.21) instead of (2.22), shows that f3 < R2S - 

Now let D := {x G R 2 : x 1 > 0, a < \x\ < (5). Then 

|^>0inD, |^ = 0in^ and L ( P) = -f^ £ < in D . 
0x1 oxi \ oxi ) oxi 
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The last inequality come from the differentiation of (1.14) and hypothesis (1.9), which implies 
that dV/dR > for r e (R-2S , R2S )- Using the monotonicity of Dirichlet eigenvalues with 
respect to the domain, this implies that /i < 0. 

Step 5 : Proof of (2.25). For any r > R, define a function fj : (r, oo) — > R by 

fj(s) := rj £ (r) exp ^a-A^s** 2 - r^^j 

where cq and p are the constants in (1.7). It follows from (2.20) and (1.7) that if s > r and 
r is sufficiently large, then V(s) — \ £ + f] 2 (s) < V(s) < c s p , so that if r is sufficiently large, 
then 

-Aft + - 2 (V -\ £ + fj 2 )fj < -Afj(s) + ^s p f) = ((-a 2 + ^)s p + + l)*^ 1 ) fj. 

Choosing a = ( - 2c °' ) ' , it follows that fj is a subsolution of (1.14) in (r, oo) if r is sufficiently 
large. For such r, noting that fj(r) = r) e (r), we can argue as in the proof of (2.21) to deduce 
that i] £ — fj is nonnegative in (r, oo). 

Then since fj(r) = r) £ (r) and fj(s) < r) £ (s) for s > r, we again use (1.7) to conclude that 

vl/2 



77' (r) > r/'(r) = — - — r^r) £ (r) > — y/2— \/v(r)r) £ (r) 

e ' e v 

for sufficiently large r. On the other hand, by choosing a = ^ in the definition of fj, 
we obtain a decreasing supersolution (still denoted fj) such that fj(r) = i] £ (r). A similar 
application of the maximum principle shows that i] £ is bounded above by (the new) fj on 
(r, 00), and in particular this implies that rf e {r) < 0. These facts combine to establish 
(2.25). □ 

We next prove 

Lemma 2.2. Assume that V satisfies (1.6) and (1.9) and the quadratic growth condition 
(1.10). Let i] £ be the positive minimizer found in Theorem 2.1. Let f £ (x) := £ e (x) /r] 2 (x), 
where £ e was defined in (1.17). Then there exists a constant C independent of e G (0, £1] 
such that 

«m)<{^' ot)+ce2/3 f:r ("a. 

In addition, for all sufficiently small e, 

HV&lloo < C (2.30) 

and 

\\fe-h\U<Ce l l\ (2.31) 
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Proof. For every s > r > R$ (where < 5 < 5o will be chosen later), we define 

ij(s) = Ve (r)e-»W* and ^ = ^ = + {K - Aq) . (2.32) 

Using (1.10), where the constant c\ is defined, and arguing as in the proof of (2.25), we find 
that f\ — r\ £ is nonnegative in (r, oo). 

We use the previous estimate and the definition of £ £ to compute 

1 P OG poo 1 

Mr) = / srf £ (s) ds < / e'^^s ds = — for r > ifc. 

The definition of / £ implies that f' £ (r) = — r — 2/ £ (r)^Jj, and from the monotonicity (2.24) 
of r) e , we infer that f' e {r) > —r in (R-s , Rs )- Thus for any R~$ < r < Rs, 

R 2 -r 2 1 

We now fix 5 = e 2 / 3 , and we conclude from (1.9) and (2.20) that (2.29) holds as long as 
r>R- So . 

For < r < R-s , we write 

T] e {rY Jr Ve\ r ) VeVJ 

From (2.22) and (1.9), we see that if < r < s < i?_ £ i/s, then 

4t4< n + ?lm2 C for sufficiently small e, (2.33) 

r) 2 (r) (1 — Ce 1 /^)^ a+(r) 

and by using the and the fact that M^-s ) < Ce 2 ^ + C5 , one easily deduces that (2.29) 
holds for r e [0,-R_ 5o ). 

Next, the definition of £ e implies that |V£ e (a;)| = \x\r] 2 (x), so that (2.30) follows from 
(2.28) and (2.21). 

For r > R_ £ i/ 3 , we see from (2.29) that |/ e (r) - / (r)| < Ce 1/3 + |/ (r)|. This is trivially 
bounded by Ce 1 / 3 if r > /?. If R < r < R £ -i/ 3 then (1.9) implies that c(R — r) < a(r) < 
C{R — r), and thus 

|/o(r)| = /o(r)<-^ f R s(R-s)ds < C(R - r) < Ce l '\ 

T — it Jr 

For < r < R_ £ i/s we write 

Mr) -Mr) = {^y J*^ s^s) ds - ^_ J^ 3 sa (s) ds^ 

, ril(R_ £ ,/z) a(R_ £ i/ 3 ) , u 

?7 £ (rJ a(r) 
= I + 11 -III 
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Using (2.33) and our earlier estimates of f £ , fo for r > R_ e i/3, we see that 

|//| < Cf e (R_ e i /3 ) < Ce 1/3 and \III\ < Cf (R_ e i /3 ) < Ce 1/3 . 
We further decompose the remaining term as 

/ 1 1 \ [ R - s l/3 1 f R -^/3 

/= —--—- / sr) 2( s )ds+—- e s(r] 2 £ (s) -a(s))ds. 

Vw) a(r)JJr a(r) Jr 

Using (2.22), it follows that 

| J| < CeV* [ R - Sl/3 s^M ds + C £ V3 / R - 1/3 d , 

Due to (2.24), < 1 if J R_ <5o < r < s < i2_ e i /3 . And if < r < i?_ <5o then ?? 2 (r) > C" 1 and 

so ^|^y < C. Thus the first integral is bounded by Ce 1 ^ 3 . The second integral is similarly 
estimated, using (1.9) in place of (2.24). □ 

Remark 2.3. In the case of a potential V for which (1.8) fails, so that for example T> has 
the form B R \ B R i, one expects that instead of being small, f e is large, namely, f £ > ce C//£ in 
the interior of B R >. This is related to the formation at very low rotations of a giant vortex in 
the interior of B R >. The arguments used to prove Lemma 2.4 show in this situation that ifV 
grows quadratically in the complement of B R , as in (1.10), then f £ is very small in H 2 \B R . 
This suggests that at low rotations there should be no vortices in R 2 \ B R , but this cannot be 
deduced from the arguments we use to prove Theorems 1.2 and 1.3. 

The last lemma in this section examines the case when V has subquadratic growth and 
f £ is also large so that in principle vortices could exist in the low density region. 

Lemma 2.4. Assume that V satisfies (1.6), (1.9) and 

there exists c 2 > and p < 2 such that V(r) < c 2 (r p + 1) for all r > R. (2.34) 
Then f e {x) —> +oo as \x\ — > oo. 

Note that with these assumptions on V, there is a sequence of functions ( a in H such that 
mi a G £ (( a ) = — oo. Physically this happens because the centrifugal force due to rotation 
is bigger than the subquadratic trapping potential. This indicates that, although one can 
prove that in this situation, f £ — > oo as r — > oo (compare Lemma 2.4), this is not expected 
to give any information about the physical behaviour of condensates. 

Proof. Let q > 2. For every r > max{l, R}, we claim that 

Ve(s) > r ]£ (r)e- u ^ q -r'>)/< ! (2.35) 

for all s > r. Where v £ r is the positive root of the polynomial z/ 2 — — £2r2g c _ 2 _ p , which 
for e small satisfy 
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with (3 = q — 1 — p/2. Indeed, the right hand side of (2.35) is a subsolution in (r, oo) of 
(1.14) while i] £ is a solution. Boths functions are going to zero at infinity and they are equal 
at s = r, so the result come arguing as in the proof of (2.21). 

We use the previous estimates and the definition of £ e to compute 

Mr) = |g = ^ fV 2 M * * f «. > £ > Cer*-*. 

and hence the result. □ 



3 Splitting the energy 

In this section we recall the proofs of (1.15) and (1.18). 

For U C R 2 , we will write E £ (w;U) etc to denote the integrals over U of the energy 
density appearing in the definition of E e {u) = E £ (u; R 2 ) , and similarly G £ (-, U), F e (-, U). 

Note that v = ujr\ £ is well defined since r\ £ > 0. Since r/ £ satisfies (1.14), we multiply it 
by r] £ (l — \v | 2 ) and integrate over a ball B r to find that 

f Br (\v\ 2 ~ + ^(V(x) + V 2 ) + i|V^| 2 ) = | fju\ 2 - V 2 ). 

Note that the Lagrange multiplier term tends to as r — > oo, since both the L 2 norms of u 
and rj £ are 1. Moreover, 

E £ (vr] £ -B r ) = J £ (r] £ ;B r ) + F £ (v,B r )+ [ h Vr] £ \ 2 (\v\ 2 - 1) + ^Vrfe • V|^| 2 

-^(1 - \v\ 2 ) 2 + ^>| 4 + ^V(x)rj 2 \v\ 2 - ^ - ^- 2 V( X W. 
We integrate by parts to obtain 

/ TiVe^Ve • V|v| 2 = - / ]-\v\ 2 Ar) £ + [ \\v\ 2 r] £ v -Vr] 

J B r I J B r 4 JdB r 2 

We use (2.25) to estimate 

| f -\ v \\u ■ V77I < - / - V 2 \v\ 2 Vv = -V(r)- 1/2 [ V \u\ 2 . 
JdBr 2 e JdBr 2 e J&B r 

Since Jr2 V|u| 2 < oo, we can easily find a sequence — > oo such that the above integral 
tends to 0. Combining the above and letting r k — > oo along this sequence, we obtain (1.15). 

The only property of V that the above argument used (implicitly) was (1.7), which will 
be used in the proof of (2.25). 
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The integration by parts that leads to (1.18) is justified in a similar fashion. One must 
estimate boundary terms of the form j dBr £v • (iv, Vi>). To do this we note that 

&-{iv,Vv) = f £ (r)r) 2 £ (ivyv) = f £ (r)(iu,Vu) < \\f e \U\u\ 2 + \Vu\ 2 ). 

We prove in (2.29) that f £ is bounded as long as V satisfies (1.10) (in fact we show that 
fe < Ce 2 ^ for large r) and since u G if 1 (R 2 ), we can again find a sequence r k — > oo such 
that the boundary terms vanish. 

Note also that the fact that f £ G L°°, or equivalently that |£ e | < Cr] 2 , implies that the 
term (, £ Jv appearing in (1.18) is integrable on R 2 for v = u/r) £ , whenever u has finite energy. 



4 Proofs of Theorems 1.2 and 1.3 



In this section we use the estimates we have already established to complete the proofs of 
our main theorems. 

Proof of Theorem 1.2. We assume that u £ minimizes E £ and that fl < C\\oge\ is such that 
(1.11) holds. 

Let x be a smooth function such that \ = 1 i n { x £ ^ : dist(x, &D) > 2| loge\~ 3 ^ 2 }, and 
with support in T>]_. We also assume that || Vxlloo < 2| loge| 3//2 . 

Let v = u £ / V£ , so that E £ (u) = G £ (rj £ ) + F £ (v) = E £ (rj £ ) + F £ (v). Thus F £ (v) < 0. We 
write 

F £ (v) = A 1 -A 2 + B 



where 



J-R? 



^v\ 2 + ^-A\v\ 2 -l) 



4e 2 



dx, 



A 2 = 2Q 



R2 



XieJvdx 



and 



B 



R 2 



(i-x) 



r, 2 £ 



\Vv\ 2 -AQf £ Jv)+^- 2 (\v\ 2 -l) 2 

loge 

r s that ft 

small £ and (recalling that | Jv\ < \\Vv | 2 ) we deduce that 

(|Vt;| 2 -4fi/ £ Jt;) > 



dx, 

It follows directly from our estimates on f £ that < f £ < C{e 2 ^ + \ loge|~ 3 / 2 ) in the support 
of 1 — x, for small enough e. Since fl < C\ loge|, it follows that flf £ < \ for all sufficiently 

2 I 



in the support of 1 — %. It follows immediately that 

-,4 



B > 



R 2 



x) 



Ve 



Vv\ 2 + ^ 2 {\v\ 2 -lf 



dx > 



(4.36) 



and hence that B = if and only if v is a constant of modulus 1 in the support of 1 — x- 
Since F £ (v) < 0, it is clear that A 1 + B < A 2 . 
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Next, define e = ^/(inf^ r] £ ), so that (in view of (2.22) and the definition of T>\) 

e<Ce\\oge\ z l\ l<i inPi. 

Then (4.36) and (2.22 imply that, 

/ ^|V^| 2 + -^(M 2 - l) 2 < (inf r l£ y\A l + 2B) < C\ log sf 2 A 2 . (4.37) 

To continue, let w = rr = w l + iw 2 . From (1.11) we see that |i>| > \ in Pi, and hence 
it is clear that w G H l (T>i), and |u>| 2 = 1. It follows that Jw = 0; we will recall a standard 
proof of this fact in a moment. Thus 

A 2 = 2Vt [ x€e(Jv - Jw) dx = 2VL ( V x (x£ £ ) • [(iv, Vu) - (iiu, Vw)] tte. 

If we write v = pe %<i> in £>!, then a calculation shows that 

(iv, Vv) = p 2 V0, (iw, Vw) = V0. 

From the latter fact we see that Jw = x (iw, Vw) = 0, as we asserted above. Also, from 
this and the fact that p > \ in V>\ we estimate 

|o 2 -l| 

\(iv,Vv) - (iw,Vw)\ = — -\pV(p\ < 2| \v\ 2 - 1| |Vul. 

p 

Using (4.37) , we deduce that 

A 2 < 2n\\V(xQ\\ooJ Vi (^\Vv\ 2 + ^(\v\ 2 -l) 2 yx 
< Cn\\V( X ^)\\oos\\oge\ 9 ^A 2 . 
One checks easily from the definitions and from (2.30) that 

||V(x6)IU < HVxIlooll^lloo + ||V6IU < C\ logef/ 2 (4.38) 

so we conclude that A 2 < Ce\ \oge\ 15 ^A 2 < \A 2 for all sufficiently small e. We know from 
(4.37) that A 2 > 0, and it follows that A 2 = 0, and hence (again appealing to (4.37)) that 
A ± — B — 0. Thus ||Vf 1 1 x,2 = ||1 — |w| 2 ||l2 = 0, and so v is a constant of modulus 1 as 
required. 

□ 

The proof of Theorem 1.3 will use the following result, which is Lemma 8 in [14]. 

Lemma 4.1. There exists a universal constant C > such that for any k G (1,2), open set 
U C R 2 andue H l (U; R 2 ), and e G ((0, 1), 

+ G E («-«/5»(l + \m„, w ) (\\<t,\\ x + 1 + / (\4,\ + l)e e (u) dx) (4.39) 
for all G C°/(U). Here e £ (u) = ||Vw| 2 + j^(\u\ 2 - l) 2 . 
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\A 2 \ < 2Qk f x ^ £ ^l-dx + S, 
Jr. 2 lose 



The lemma as stated in [14] does not explicitly specify the exponent (/c — 1)/50 appearing 
on the right-hand side of (4.39). By inspection of the proof, however, one sees that this 
exponent can be taken to have the form |a, where a — (k — l)/12/c as in Theorem 2.1 of 
[15]. 

Proof of Theorem 1.3. We continue to use notation from the proof of Theorem 1.2, such as 
Al, A 2) B, e, and so on. 

We first invoke the lemma, with e in place of e and x £,e m place of 0, and with k > 1 to 
be chosen. This yields 

logs | 

where S denotes the error terms in (4.39). We note that for all sufficiently small e > 0, the 
error term satisfies the bound 8 < Ce^(l + \A 2 \), for f3 — (k — 1)/100, for all sufficiently 
small e. This is a consequence of (4.37) and the estimates 

\\x &\\w^ < C\ loge| 3/2 , ||x < c. 

These in turn follow from (4.38) together with (2.30) . 

Now the choice of e implies that eg(v) < ^|Vv| 2 + ^(|t>| 2 — l) 2 inV 1 , and recalling that 
Ce = feVei we obtain 

(l-Ce?)\A 2 \ < 2^M|/ x (f|V,| 2 + ^(|,| 2 -l) 2 ) + ^ 

I logs I 

We know from (2.31) that H/J^ < (l + Ce^H/olloo < (l + Ce^H/olU, and from the choice 
of e, for any > there exists e > such that | log e| > (| loge| — log | loge|)(l + Ke 13 ) if 
< e < e . Thus 



\A 2 1 < n ( ^1122- ) «A X + 

\\\oge\-\og\\oge\J 

for all sufficiently small e. Assume that Q < 2 y^ (| log e| — (ci + 1) log | log £ | ) , for Ci to be 
chosen below. Then 



log 


loge 




I logs: 




log | 


logs] 



y Al + CeP < ^-c ^I^^ kA. + Cs^. (4.40) 



\A 2 \ < \i- Cl 

We now take k := 1 + c x lo ^°* £ 1 , so that = (k - 1)/100 = jfc lo ^y 1 . Recalling that 
A 1 + B < A 2 and that B > 0, clearly A x < A 2 , so we deduce that 

ci( ^zfl y Al <Cee = c\i oge\-^™. 
\loge\ 

If c\ = 400 then we conclude that A\ < C\\oge\~ 2 . 
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Then (4.40) implies that A 2 < C\ \oge\ 2 , and it follows that B < C\ log e| 2 . In view of 
(4.37), this implies that 

/ |Vtf + -^(M 2 -1) 2 <C|log £ |- 2 . (4.41) 
JVi 4e z 

The estimate ||Vi>||oo < 7 (see (2.23)) and (4.41) are easily seen to imply that 

\v\ > 1 — C| logel" 1 in V 1 (4.42) 



for all sufficiently small e. Thus Q is subcritical for small enough e. 



□ 
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